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MAGNETIC RESONANCE IMAGING SEQUENCE DESIGNER 

BACKGROUND OF THE INVENTION 
Technical Field of the Invention 

The present invention relates to the field of magnetic resonance imaging (MRI) and, 
in particular, to a system and method for aiding the efficient design of pulse sequences. 

Description of Related Art 

Until the development of MRI and Nuclear Magnetic Resonance (NMR) 
technology by Dr. Raymond V. Damadian in the 1970*s, diagnostic imaging of internal 
physiology was limited to techniques which provide limited soft tissue contrast. For 
example, as is well understood in the imaging art, computed tomography (CT) techniques 
depend on tissue density, e.g., soft tissue compared to bone, and usage of contrast media, 
e.g., barium, both affecting x-ray attenuation and detection. Although CT, at present, 
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reveals better bone detail, MRI is far superior for most other soft tissues, illuminating the 
internal networks and pathways to physicians without the known deleterious effects of x- 
rays. 

Although a full description of how MRI works is not necessary to the understanding 
5 of the subject matter of the present invention, a brief illustration of the physical principles 
involved is set forth below. In short, MRI is a diagnostic method for providing detailed 
specimen images through manipulation of atomic nuclei, specifically hydrogen, within a 
specimen tissue. A fundamental property of individual nuclear particles is that individual 
particles spin or rotate about their own respective axes. As is understood in physics, a 

1 0 spinning charged particle produces a magnetic moment directed along that particle's axis of 
rotation. These spinning nuclei and their resulting moments are randomly oriented in the 
absence of any external magnetic fields. However, by applying a magnetic field, the 
rotating nuclei essentially align their axes either in parallel or in opposition to the magnetic 
field. Those nuclei aligned in opposition to the magnetic field have a higher energy than 

15 those nuclei that are aligned in parallel with the field. A small majority of nuclei will be 
aligned in the lower energy state, i.e., in parallel, than opposed to the same field, usually 
only measuring in parts per million for the excess. By the addition of energy, e.g., by 
application of radio frequency (RF) energy, to these lower energy state excess nuclei, these 
nuclei can be transitioned to align themselves antiparallel or in opposition to the magnetic 
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field. As is understood in the art, it is these few reaHgned nuclei that ultimately provide the 
information used to generate an MRI image. 

While the respective nuclei are generally aligned with the applied magnetic field, it 
should be understood that this alignment is not precisely with a plane parallel to the axis of 
5 the magnetic field. Instead, the nuclear moments align at a slight angle from the axis of the 
magnetic field and precess about this axis. This frequency of precession, along with the 
magnetic moment caused by the alignment of the nuclei, comprise the phenomenon on 
which imaging by magnetic resonance is based. 

The frequency of this atomic or nucleic precession, also referred to as the Larmor 

10 frequency, is a function of the specific nucleus and the strength of the external magnetic 
field. The nuclei will absorb energy and induce a signal in adjacent RF receptor coils only 
at the particle's Larmor frequency - an event referred to as "resonance." In other words, by 
applying energy to the specimen at the Larmor frequency, the net magnetic moment of the 
excess nuclei may be reversed, or deflected, to the opposite or antiparallel direction by 

15 causing these parallel state particles to elevate to the higher energy state. The 
radiofrequency energy pulses applied are referred to as "excitation pulses." The duration of 
the RF pulse specifies the duration of the nuclear moment deflection. When the excitation 
pulse is removed, the nuclei will then begin to lose energy, causing the net magnetic 
moment to return to its original, lower energy state orientation, and the energies emitted 

20 during this transmission are used to create the image of the specimen. 
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Present day MRI devices generally scan only hydrogen atoms. The hydrogen atom 
is most attractive for scanning since it comprises the largest atomic percentage within the 
human body and provides the largest magnetic resonance (MR) signal respective to other 
elements present in human organs. As described hereinabove, every nuclear particle spins 
5 about its axis and the individual properties of the spin are defined by the specific nuclear 
particle in question, e.g., hydrogen, creating a magnetic moment v^ith a defined magnitude 
and direction. The Magnetic Resonance (MR) signal itself is a complex function 
dependent upon the concentration of the deflected hydrogen atoms, spin-lattice relaxation 
time (Tl), spin-spin relaxation time (T2), motion within the sample and other factors as is 

10 understood in the art. 

Another component of the MR signal is, of course, the particular series of RF and 
magnetic field gradient pulses employed in the form of pulse sequences. Varying the pulse 
sequences can produce considerable image differences, such as Tl emphasis (Tl- 
weighted), T2 emphasis (T2-weighted), proton density emphasis or combinations thereof. 

15 Common sequences include Gradient Echo (GE), Spin Echo (SE), Inversion Recovery 
(IR), Double Spin Echo, 3-dimensional Gradient Echo (3DGE), 3-dimensional Spin Echo 
(3DSE), Fast Spin Echo (FSE), Partial Saturation (PS) and others. It is understood that 
these sequences are illustrative only and the present invention is in no way limited to 
application of only these specific sequences. Since one sequence image type may not 
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optimally illustrate an area of consideration, multiple images using varying sequences of 
pulses may be required to fully analyze the area, as is understood in the art. 

At present, conventional MRI systems offer fairly primitive interfaces for the 
design of the aforementioned pulse sequences. In particular, present MRI systems are ill- 
suited for sequence designers who must input and modify customized pulse sequences. 
Furthermore, this input is generally made by coding the sequence in a programming 
language, e.g., C, and is further complicated in that the coded sequence format must be 
tailored for each individual machine, thus necessitating that the sequence designer must be 
skilled in the programming arts along with the MRI technologies or alternatively requiring 
an MRI sequence designer to work in coordination with a programmer. Consequently, 
conventional systems generally lack real-time communications with the MRI scanner since 
each sequence must first be coded and compiled prior to being loaded on the system. 

Accordingly, a first object of the present invention is to provide an improved MRI 
apparatus for more efficient creation and development of MRI pulse sequences. 

It is a second object of the present invention to provide a graphical user-interface 
for performing the mathematical calculations related to MRI pulse sequence design, 
thereby providing an automatic graphical response of the interface to user manipulation, 
facilitating interaction between the user and the interface. 

It is a third object of the present invention to provide a graphical user-interface for 
intermediating between the sequence designer and the MRI hardware such that the designer 
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can directly view the details of the entire pulse sequence but can also access and modify 
the sequences directly through a mouse or keyboard. 

It is a fourth object of the present invention to provide a real-time interface, or 
front-end, between the graphical user-interface and the MRI system hardware that enables 
real-time communication and interaction between the sequence creator and the MRI system 
hardware. 

It is a fifth object of the present invention to provide real-time communication and 
interaction between the sequence creator and the MRI system hardware, enabling data 
acquisition and graphical display of the RF shapes, gradient waveforms and MRI signals 
received inside the magnetic field and providing analysis of this information in real-time. 

It is a sixth object of the present invention to provide a real-time communication 
and interaction between the sequence creator and the MRI system hardware, thereby 
enabling dynamic manipulation of the details of a MRI pulse sequence accessed through 
the graphical user-interface. 

It is a seventh object of the present invention to provide a real-time communication 
and interaction between the sequence creator and the MRI system hardware that enables 
detection of dynamic deficiencies of the MRI system through feedback information and 
possible compensation for the deficiencies through sequence manipulation. 

It is an eighth object of the present invention to provide a foundation for 
development of automated calibration of imaging sequences of an MRI system. 



PATENT APPLICATION 
Docket No. 459713-1 



SUMMARY OF THE INVENTION 

The present invention is directed to a system, apparatus and method for facilitating 
magnetic resonance imaging (MRI) pulse sequence generation and modification and for 
5 real-time sequence input modification for use in conjunction with magnetic resonance 
imaging equipment. A graphical user interface is provided through a display coupled to a 
digital computer operating as the primary control system for a magnetic resonance imaging 
scanner and associated hardware. Through the graphical user interface, an operator, either 
at the site of the MRI unit or at a remote location, may choose or design sequences of 
10 radiofrequency pulses, gradient waveforms and other input parameters for the magnetic 
resonance imaging apparatus. Real-time information is also communicated to the operator 
through the graphical user interface allowing for real-time manipulation of the magnetic 
resonance imaging inputs and for displaying the magnetic resonance response thereto. 

15 BRIEF DESCRIPTION OF THE DRAWINGS 

A more complete understanding of the system, method and apparatus of the present 
invention may be obtained by reference to the following Detailed Description when taken 
in conjunction with the accompanying Drawings wherein: 

FIGURE 1 is a simplified diagram of a magnetic resonance (MRI) imaging system 
20 upon which the principles of the present invention are performed; 
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FIGURE 2 illustrates an exemplary MRI scan controller interfacing the control 
system with the MRI hardware on which the present invention can be applied; 

FIGURE 3 illustrates an exemplary graphical user interface from which a user can 
initiate a sequence design according to a preferred embodiment of the present invention; 

FIGURE 4 illustrates an exemplary settings dialog box for accepting user defined 
resolution parameters; 

FIGURE 5 illustrates an exemplary sequence-type dialog box for allowing a user to 
define various parameters associated with a given type of pulse sequence; 

FIGURE 6 illustrates an exemplary sequence-tailor dialog box providing various 
controls for user selection; 

FIGURE 7 illustrates an exemplary sequence display for displaying various 
sequence plots to the user and allowing user interaction therewith; 

FIGURES 8 A through 8C provide additional illustration of an exemplary sequence 
display having additional control elements invoked for assisting a user in timing analysis of 
the various plots; 

FIGURES 9A through 9C illustrate various control dialog boxes available through 
invocation of a shape editor command; 

FIGURES lOA, lOB, and IOC illustrate an imaging parameters dialog box for 
displaying and accepting modifications to various imaging parameters; 

FIGURE 1 1 illustrates a scan gradient wave option box and displayed parameters; 
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FIGURE 12 is an exemplary graphical representation of a variable window used in 
accordance with the Gradient Wave Scan dialog box of FIGURE 11; 

FIGURES 13 A and 13B illustrate plot editing commands available for modifying 
and design of a pulse sequence; 
5 FIGURES 14A and 14B illustrate an exemplary embodiment of integral condition 

calculations according to a preferred embodiment of the present invention; 

FIGURES 15 A and 15B provide further illustration of the automated integral 
conditioning of the present invention; 

FIGURES 16A through 16D provide illustration of acceleration compensation in a 
10 slice select (SS) gradient, velocity compensation in signal recognition (RO) and automatic 
integral conditioning thereof; 

FIGURE 17 illustrates a secondary window invoked at real-time interfacing 
initiation; 

FIGURES 18A and 1 8B illustrate an exemplary tuning display and illustrative plots 
15 displayed therein; and 

FIGURE 19 illustrates an exemplary scan parameters dialog box according to a 
preferred embodiment of the present invention. 
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DETAILED DESCRIPTION OF PRESENTLY PREFERRED EXEMPLARY 
EMBODIMENTS 

The numerous innovative teachings of the present application will be described 
5 with particular reference to the presently preferred exemplary embodiments. However, it 
should be understood that this class of embodiments provides only a few examples of the 
many advantageous uses of the innovative teachings herein. In general, statements made in 
the specification of the present application do not necessarily delimit any of the various 
claimed inventions. Moreover, some statements may apply to some inventive features but 
10 not to others. 

As discussed hereinabove, current magnetic resonance (MR) imaging apparatuses 
may utilize a number of different pulse sequences. More particularly, MR images are 
obtained by using an appropriate sequence of specific RF pulses, signal (echo)-gathering 
times (TE) and sequence repetition times (TR). For example, dependent on the desired 

15 image emphasis, e.g., Tl, T2, or proton density, specific sequence types can produce 
dramatically different imaging results. As discussed, examples of common pulse 
sequences include Gradient Echo (GE), Spin Echo (SE), and Inversion Recovery (IR), as 
well as Double Spin Echo (DSE), 3-D Gradient Echo (3DGE), 3-D Spin Echo (3DSE), and 
Fast Spin Echo (FSE). Again, it is understood that the present invention is not limited by 

20 the sequences listed above. The present invention greatly simplifies the manipulation of 
MR imaging parameters by allowing for more efficient sequence design and parameter 
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tailoring via a digital computer, as discussed in more detail hereinbelow. The present 
invention allows a designer to generate and modify a variety of sequences quickly and 
efficiently through a graphical user interface coupled to a digital computer, where the 
digital computer is itself may be coupled to the magnetic resonance imaging equipment. 
Additionally, the system of the present invention allows for real-time communication with 
the MR scanner providing real-time viewing of MRI signals, gradient waveforms and RF 
(pulse) shapes. Furthermore, sequences can be calibrated automatically and dynamically in 
response to parameter input from the designer. 



The present invention may be more readily understood with reference to FIGURE 1 
10 in which an exemplary magnetic resonance imaging (M^I) unit 150 as implemented in a 
preferred embodiment of the present invention is generally depicted. A specimen aperture 
155 is centrally located within the bulk of th^MRI unit 150 scanning apparatus. 
Surrounding the specimen aperture 155 are gradient coils 180, transmitter coils 190 and 
receiver coils 195. The gradient coils 180 provide the gradient waveform for facilitating 
15 selective excitation, echo formation and/localization, among other purposes, as is 
understood in the art. Generally, the gpdient coils are responsible for generating the 
gradient field in corresponding axes iff mutually orthogonal directions, e.g., designated 
according to the Cartesian x, y and z Coordinates, and thus are driven independent from one 
another. The transmitter coil 19(9 is responsible for delivering the RF pulses into the 
20 aperture 155 at the Larmor frequency for providing excitation of the protons in the 
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specimen, e.g., a human body. The receiver coil 195 is responsible for acquiring the signal 
generated from the excited protons during the proton relaxation period, as is understood. 
Driving these coils and interfacing the MRTiinit 150 with a computer 110 is a scan 
controller 130. The scan controller 130is the hardware that provides the real-time data 
delivery to and from the various e^ipment, e.g., digital to analog converters, as described 
more fully hereinbelow wi^reference to FIGURE 2. An iron yoke 170 completes the 
circuit by coupling th^oles of the magnetic material of the plates 200. It is understood 
that permanent magnets, as well as superconducting magnets, may be suitably interchanged 
for the resistwe electromagnet. 



□ 10 Magnetic material provides a primary static, i.e., uniform and constant with respect 

in 

H to time, magnetic field for surrounding the specimen to be imaged. Gradient fields 

O corresponding to the Cartesian coordinates are used for coding position information with 

respect to the MRI echoes. Therefore, three separate gradient coils 180 are required, each 
coil being independently driven. A digital computer 110 equipped with standard input 
15 devices, e.g., a keyboard 1 15, a mouse or pointer device 140, and output devices, e.g., a 
display 120, facilitate the sequence designer's interaction with the scan controller 130 and 
thus the overall MRI unit 150. An MRI pulse sequence containing all the information, e.g., 
three gradient waveforms, RF pulse shape definitions, signal acquisition timing data, etc., 
for generating an MRI signal, i.e. an echo, can be designed, modified or stored within 
20 computer 1 10 for application to the system according to the present invention. 
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With reference now to FIGURE 2, a more detailed illust^ion of the scan controller 
210 interfacing between computer 200 and MRI unit 220 is ilmstrated. Scan controller 210 
is in communication with computer 200 over data hnk:s^02 and 204 for respectively 
sending and receiving digital data thereover, although a smgle bi-directional data link could 
5 be substituted therefor. Timing of all interactions, including frequency and phase settings 
for the respective transmitter 190 and receiver coil/95, with MRI unit 150 is coordinated 
by synthesizers 230 integrated within scan controller 210. In a preferred embodiment, 
synthesizers 230 include two synthesizers which interface with two synthesizers DACs 23 1 

/ 

and 232. The scan controller 210 interfaces with three gradient DACs 240, 242 and 244, 

DACs 250 and 252. To control and operate an advanced MRI unit as in the preferred 
embodiment, at least seven DACs are needed to transfer digital input between the scan 
controller 210 and the MRI uni^ However, the preferred embodiment is capable of 

/ 

numerous modifications and rearrangements which would require more or less than seven 

/ 

1 5 DACs. Gradient DACs 240 are supplied with digital gradient information describing the x, 

/ 

y and z gradient fields over respective gradient DAC input lines 241, 243 and 245, as 
defined by the MRI pulse sequence received by scan controller 210 from computer 200. 
Gradient DACs 240, 242 and 244 then convert the digital gradient data to corresponding 
analog gradients wWch are transmitted to respective gradient coils 1 80 over gradient DAC 
20 output lines 246,^47 and 248. Thus, gradient DAC 240 is responsible for driving one 
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gradient coil, e.g., gradient coil ISOX^ofthe set of gradient coils 180. Likewise, each of 
gradient DACs 242 and^44are responsible for driving a corresponding one of the 
remaining gradipifcoils 180, e.g., respective gradient coils 180 Y and 180 Z. 



RF shaper DACs 250 and 252 are responsible for converting the RF shape data 
5 from the digital domain as defined in the MRI pulse sequence received by the scan 
controller 210 from computer 200 and /modulating the representative RF shapes 
accordingly. The RF shapes are received/and modulated by transmitter coil 190, generally 
at the Larmor frequency, to the subject specimen being analyzed in aperture 155. The 
frequency and phase of these modulations are controlled by synthesizers 230 over control 

1 0 lines 260 and 26 1 . The echo resi/lting from the applied gradient waveforms and RF pulses 
is acquired by the receiver/coil 195 during the relaxation periods and accordingly 
transmitted to the scan ccmtroller 210. Proper acquisition of the echo is facilitated by 
frequency and phase settings applied to the receiver channel by synthesizers 230. These 
frequency and phas/settings are supplied to the scan controller by the MRI pulse sequence 

15 data received frcyh computer 200. Thus, the MRI pulse sequence provided by computer 
200, under cojnmand of an operator thereof, directs the operations for echo generation and 
acquisition) 



20 



Digital computer 110 includes a processor, e.g., a microprocessor from the family 
of Pentium processors manufactured by the Intel corporation for directing and 
performing operations and receiving and executing input from a user, e.g., from the 
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keyboard 1 15 or pointer device 140. Digital computer 1 10 also contains a bank of random 
access memory (RAM) for storing and executing commands therefrom, and a long-term 
storage media, e.g., magnetic disk, for storing executable instructions that are retrievable 
and loadable into RAM. In a preferred embodiment, digital computer 110 has a Microsoft 
5 Windows™ operating system for coordinating and executing instructions and programs, 
coordinating communications to peripheral hardware, scheduling tasks and allocating 
hardware. The present invention preferably includes a Windows executable program 
stored in long-term storage media and executable from RAM, although other platforms are 
not precluded. 

10 The present invention allows for efficient creation and customization of generic 

pulse sequences through a primary design interface 300, i.e., a^^^hical user interface 
presented in the form of a window 305, presented on dispjay i20 and generally depicted in 
FIGURE 3. The design interface consists of a numfe^ of user selectable menu editor items, 
e.g., New Sequence 310, Sequence Tailor 320, and Delta Size 330, as well as the common 
15 Windows™ menu editor items, e.g.^^Hle 340, Edit 350 and View 360, all contained within 
the mainframe window 305^,/Various toolbar editor items, e.g., Open File 370 and Save 
380, may be locatedJfefQow the menu editor items. Selection of a given menu choice 
initiates generprion of a drop-down box, e.g., drop-down box 390 resulting from user 
selectioj^f the new sequence menu editor item 310. Each of the menu editor items will 
20 g^rally have a drop-down box associated therewith, each of which generally provide the 
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user with additional selections associate^ith the selected menu item. A specific dialog 
box may be activated by selectioiLOT a drop-down box item. For example, a Nul Sequence 
, drop-down box item 396, wfien selected, activates a specific dialog box (not shown) that 
allows the user to cmme a sequence from scratch. The user has unrestricted freedom to 
build a sequence from existing types, while providing the advantages of copy, 
modifi^^don, and other functions provided by the user interface. 

In FIGURE 4 is illustrated an exemplary Settings dialog box 400 displayed within 
the mainframe window 410 for accepting resolution data for a new sequence. Settings 
dialog box 400 preferably is activated by selection of any option within drop-down menu 
390. The Setting dialog box 400 preferably includes an RF resolution box 402, a gradient 
resolution box 404, and a graphic resolution box 406. A user may enter a desired 
numerical value, as depicted in this illustrative example, or accept the system defaults, e.g., 
4 and 16 microseconds for the respective RF and gradient settings. The graphic resolution 
box 406 allows the user to have an option between speed and clarity of the graphic 
drawing. The user may choose the default graphic resolution, or may opt to enter a higher 
or lower resolution. When entering a higher resolution, screen repainting requires more 
time than the default resolution. The lack of speed becomes an issue when performing a 
sequence with long train echoes, such as FSE and EPI with 127 or 256 echoes. The 
displayed settings are accepted by selecting the OK button 408. 
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A sequence parameters dialog box 500, as illustrated in FIGURE 5, is displayed in 
the event that the resolution settings are accepted in the previously displayed Setting dialog 
box 400. The sequence parameters dialog box 500 displays various user definable attributes 
or parameters associated with a selected generic pulse sequence. As illustrated in FIGURE 
5 5, the sequence parameters dialog box 500 preferably includes default values, such as a 
sequence name box 502 and a TE field box 504 indicating 85 ms, which are prestored with 
the associated generic waveform characteristics within computer 110. A number of 
attribute fields are also included in the sequence parameters dialog box 500, a group of 
gradient motion compensation field values 506, various RF pulse characteristics field 

10 values 508, and data acquisition field values 510. A user designing a sequence can, of 
course, accept the default values by selecting the confirmation command, e.g., clicking on 
an OK button. Doing so results in a logical association of the depicted parameters with 
respect to the selected generic pulse sequence. The user also has the option of canceling 
the sequence design by selecting a cancel button. Advantageously, the user also has the 

15 option to manually modify the displayed sequence or acquisition parameters by manually 
entering a desired value in any one of the sequence parameter boxes. 

FIGURE 6 illustrates a sequence tailor dialog box 600 that is preferably displayed 
after acceptance of the sequence parameters from sequence parameters dialog box 500. 
The sequence tailor dialog box 600 comprises four main control features in which the user 

20 may interact: controls section 602, shape editor 604, block editor 606, and time scaling 
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608. Each of these sections have at least one object through which the user may supply 
input. These objects may include the calibration check box 610, delta on check box 612, 
the update lock checkbox 614, and the shape editor radio button 61 6, or any other graphical 
object useful for obtaining user input. 
5 The various controls activated through the sequence tailor dialog box 600 are 

preferably available for user interaction therewith during which the current sequence 
design is displayed to the user according to an exemplary sequence display 700 as 
illustrated in FIGURE 7. The sequence displayed corresponds to a sequence selection 
made as described hereinabove. The sequence name as entered in the sequence name field 
10 502, is provided on the sequence display 700 title bar 702. The graphical display of a 
sequence within the sequence display 700 is preferably divided into four general sections. 
The topmost section of sequence display 700 provides a graphical display of the RF 
characteristics of the sequence as currently designed and is designated as the RF display 
704. 

15 A second portion 706 of the sequence display 700 of FIGURE 7, displays a slice 

select (SS) gradient graph 708, representing the gradient used for a particular slice selective 
excitation. For example, if an axial (transverse) specimen image is desired, the slice select 
gradient 708 would be termed "Gz", meaning that the external field is aligned with the z- 
axis. It is understood that slice select gradient may be taken on any of the Cartesian axes, 

20 such as Gx or Gy, depending on the desired slice orientation. 
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The third portion, RO section 710, of the sequence display 700 of FIGURE 7 
represents the dephasing or signal acquisition (RO) graph 712. The RO section 710 
essentially provides signal location information for a given echo and is typically calculated 
by a fast Fourier transform algorithm. Final image generation may be performed locally by 
5 the digital computer 1 10 or may be offloaded to high speed numerically intensive systems. 
'■^ The lowermost portion 7 1 4 of the sequence display depicts the phase encoding (PE) 

1% graph 716. 

fii 

i3 Four horizontal lines 718, 720, 722, and 724 respectively indicate the zero- 

Q 

lij amplitude of the respective RF section 704, SS section 706, RO section 710, and PE 

Q 10 section 714. The time coordinate commonly shared among each of the plots is represented 

according to standard convention along the horizontal and originating from the leftmost 
side of the sequence display 700. Preferably, the four sections of the sequence display 700 
are automatically scaled according to calculated maximum and minimum amplitudes of the 
w^aveforms displayed therein. These calculations are performed upon confirmation of the 
1 5 original setting by selection of the OK button in settings dialog box 400. The conventions 
preferably defining these scaling calculations are: 




RF 580 - a maximum positjj^d^amplitude of a RF shape is of 100 scaling units; 
SS 704 - the gradienUflateau corresponding to the slice selective RF pulse is of 100 
scaling units; 
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5 

ia 10 



RO 706 - the gradient plateau corfesponding to the data acquisition window is of 
100 scaling units; and X 



\f PE 710 - the absoluj^^^aximum amplitude among all the PE plateaus, which are 
stepped during the sc^is of 100 scaling units. 

Detailed information regarding a particular point of a given plot may be obtained 
through user interaction with the user interface 305 preferably through directions of the 
mouse 140. An exemplary tool tip box 726 is displayed when the mouse pointer is 
positioned over a given point, or node of a displayed plot. As illustrated, the tool tip box 
726 provides detailed numerical data representative of the subject node. The exemplary 
tool tip box 726 indicates that the selected node is the tenth node (starting from node zero) 
along the RO gradient and represents a timing of 43.52 milliseconds along the plot. 
Furthermore, the RO amplitude is also provided (0.000) as well as the respective time 
differences between the previous node (1.024 ms) and the next node (2.512 ms) of the 
associated plot. Data related to the other plots can be obtained by simply directing the 
mouse pointer to a displayed node on any of the displayed plots. 

As previously mentioned with respect to FIGURE 6, sequence tailor dialog box 600 
is provided with various controls such as Block Editor Single 618 and Global 620 radio 
buttons and time scaling Xscaling 622 and Normalize 624 radio buttons, and preferably is 
displayed concurrently with user interface 305. 
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These various controls can be utilized for further plot enhancements, as illustrated 
in FIGURES 8A-8C. It should be understood to one skilled in the art that the 
aforementioned Xscaling is for the purpose of inserting new waveforms at the beginning, 
middle, or end of an existing waveform. Various vertical lines 802, 804, 806, 808, 810, 
5 812, 814, 816, and 818 are displayed in response to user selection of the Xscaling radio 
button 622. These vertical lines pass through each of the displayed plots and assist the user 
in analysis of the timing relations among the RF pulse, data acquisition and gradient 
waveforms. An example of such timing relation is indicated in FIGURE 8B. The duration 
of the slice plateau should be approximately equal to or slightly larger than the 

10 corresponding RF shape duration. In this illustrative example, a user can quickly make an 
assessment of this relation by selecting the Xscaling radio button 622 and noting the 
relation of the RF shape duration spanning vertical lines 802 to 806. In this instance, the 
SS plateau duration is slightly larger than the RF pulse duration, indicating a proper 
relationship. Various other timing relations can easily be confirmed, such as verification 

15 that a PE gradient occurs only during periods with no RF pulse or data acquisition, by 
viewing the exemplary plots as graphically displayed in the preferred user interface 305. 
FIGURE 8c is illustrative of the vertical scaling dependency on the subject sequence. 

Further controls are provided by selection of the shape editor radio button 616 in 
the sequence tailor dialog box 600 as depicted in FIGURE 6. When this control is selected, 

20 a shape editor dialog box 900 is displayed as illustrated in FIGURE 9A. This dialog box 
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allows the user to either modify the RF shape parameters by selection of the modify 
parameters button 902, or modify the RF shape itself by selection of the select types button 
904. 

Selection of the modify parameters button 902 generates a shape modification 
5 dialog box 910, as illustrated in FIGURE 9B. The shape modification dialog box 910 
displays the current settings of the parameters of the selected RF shape. A user can modify 
the parameters associated with the current RF shape by changing the parameters to the 
desired value, entered by the user, of any of the displayed parameters such as amplitude of 
sine 912, cycles of sine 914, exponential of sine 916, amplitude of cosine 918, cycles of 
10 cosine 920, exponential of cosine 922 and asymmetric factor 924. Upon acceptance of the 
modified parameters, the RF plot, SS plot 708, RO plot 712, and the PE plot 716 of 
FIGURE 7 are redrawn for display accordingly. 

If the user selects the select types buttpn 904 in the shape editor dialog box 900, as 
shown in FIGURE 9A, a shape selection/malog box 950 is displayed. A shape drop down 
,1 5 box 952 within the shape selectiondmlog base 950 displays the current selected shape. By 
selecting the drop down arrow^54, the shape drop down box 952 expands to provide the 
user with a selection oL^vailable RF shapes, such as Gausian, Gasico, Hermite, etc. 
Selection of a RF sb^e different than the current RF shape will result in the RF plot, SS 
plot 708, RO p)^ 712, and the PE plot 716, being redrawn for display upon return to user 
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terface 305. Preferablyy^restored numerical shape files may be accessed through 
selection of the shapes available via button 956. 



When a sequence is created and modified, it can be saved with the original user- 
assigned name by selection of the Save item (not shown) provided by selection of the File 
5 menu editor item 340, as illustrated and described in connection with FIGURE 3. 
Alternatively, a new name can be assigned by choosing the Save As item from the File 
menu editor item 340. In either case, the saved sequence is digitally stored in, preferably, a 
plurality of files, each of which defines a particular aspect of the saved sequence. These 
files are preferably stored within a storage device of computer 110, as illustrated in 

10 FIGURE 1. In a preferred embodiment, there are eleven files created and saved when a 
user chooses to save a given sequence. These files include seven binary files for input to 
the various DACs described hereinabove, as well as a documentation file, e.g., a MS Word 
document (.doc), a source file such as a C++ file (.cpp) and a graphical file for sequence 
plotting such as a .gwa file. The seven binary files for directing the DAC operations can be 

15 combined into a single file to reduce the number of stored files. The documentation file 
preferably provides all the numerical values describing the sequence including a report on 
various critical values related to the sequence and the overall system, e.g., maximum rising 
time and electrical currents of the associated gradient. 

After a sequence is designed and stored, an MRI scan is available to be performed. 

20 This is initiated through selection of the ScanSettings menu editor item 1 000, as illustrated 
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in FIGURE lOA. The ScanSettings menu editor item 1000 is included in the main menu of 
the primary design interface 300 but only the ScanSettings 1000 menu editor item and its 
associated drop down menus 1002, 1020, and 1040, and a Setting Combo Imaging 
Parameters dialog box 1050 A are illustrated in FIGURE lOA. When the ScanSettings 
menu editor item 1000 is selected, the ScanSettings drop down box 1002 is displayed and 
provides the user with a selection of scan parameters to perform a scan. 

In the illustrative example, the MRI Scans option 1004 is selected and causes the 
MRI scan type drop down menu 1020 to be displayed. Eight different types of scans are 
illustratively available for selection, e.g., 2D-Scan 1022, 2D- Variable TR Scan 1024, 3D- 
scan 1026, FSE 2D-scan 1028, FSE 3D-scan 1030, Multiple 2D-Scans 1032, Multiple 3D- 
Scans 1034, and Combo Scan 1036. It should be understood, however, that the number 
and type of available scans is not limited to the illustrative examples depicted. 

If a scan type selected does not have an associated sequence type already loaded as 
aforedescribed, an error message is generated and displayed to the user indicating the 
absence of the desired sequence type. Additionally, the error message is generated to 
remind the user that the desired scan did not match the type of sequence. A simple 
example of the image parameter setting is illustrated in FIGURE IOC. A Setting Imaging 
Parameters dialog box 1080, as illustrated in FIGURE IOC, is displayed upon selection of 
MRI Scans type having an associated sequence type already currently loaded. The Setting 
Imaging Parameters dialog box 1080 also displays numerous other parameters which are 
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listed above with reference to FIGURE lOA. The default settings are displayed consistent 
with the sequence already displayed in the user interface 305. In the illustrative example of 
FIGURE lOA, the Combo Scan option 1036 is selected, which activates the display of a 
combination type drop down menu 1040. The user may want to perform the combo scan 
5 1036 when performing a coronary MRA or MRI guided surgery. In the preferred 
embodiment, a user may choose between a StaticCombo scan 1042 and a DynamicCombo 
scan 1044. When the StaticCombo scan 1042 is selected, as illustrated in FIGURE lOA, 
the Setting Combo Imaging Parameters dialog box 1050A is displayed. In this illustrative 
example, the first sequence denoted in a first sequence name box 1052 is used for an 

10 extremely fact detective scan, whereas the second sequence denoted in a second sequence 
name box 1054 is used for actual imaging. The name of the combination of sequences is 
preferably provided in a sequence name box 1056. A default first sequence echo gathering 
time parameter 1058, and a default second sequence echo gathering time parameter 1060 
are also displayed. Various other imaging parameters are preferably provided, e.g., number 

15 of slices 1 062, slice thickness 1 064, sequence repetitions, number of phase encoding levels 
1066, Discrete Fourier Transform (DFT) size 1068, and 180 RF polarity flipping. These 
displayed parameters are exemplary only and can be tailored according to the software 
designer's preference. The numerical values associated with the various displayed 
parameters are modifiable by the user by simply entering different values from the defaults 

20 and accepting the entered values by selection of the OK button. Upon acceptance of the 
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default values, or values substituted therefor, two files (.va and .exam) are preferably 
generated and saved accordingly. These files provide information for proper setting of 
sequences, scan and image reconstruction required for the various imaging components. 

A Setting Combo Imaging Parameters dialog box 1050B, as illustrated in FIGURE 
lOB, is displayed when the user selects to perform a DynamicCombo scan 1044. The 
trigger sequence used is denoted in a trigger sequence name box 1070, and the imaging 
sequence used is denoted in a imaging sequence name box 1072. Various other 
parameters, such as those listed above with reference to FIGURE lOA, are also displayed 
in the Setting Combo Imaging Parameters dialog box 1050B. 



10 FIGURE 11 illustrates a preferred embodiment o^e present invention when a 

Scan a Gradient Wave option 1006 is selected frojrithe ScanSettings menu editor item 
1000. Upon selection of the Scan a GradienLWave option 1006, an interference dialog box 
1 100 is displayed. Selecting the OKJalitton on the interference dialog box 1 100 causes a 



3^ 



Setting of Gradient Wave Scan^alog box 1120 to be displayed. The Setting of Gradient 
15 Wave Scan dialog box W70 allows a user to enter a particular window location and size 
with which to copa^are the input and output retrieved from the magnetic field. The user 
may entery^ues in a variety of parameter boxes, e.g. a sequence name box 1122, a 
sequ^e repetition times box 1 124, a number of data acquisition box 1 126, and a sample 
terval box 1128. 
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FIGURE 12 illustrates an exemplary graphical representation of a variable window 
1200 which can be selected in the Setting of Gradient Wave Scan dialog box 1120 in 
FIGURE 1 1 . In this example, the user selected the SS gradient with a sample interval 1 128 
of 16 microseconds. One skilled in the art will recognize that the settings pertaining to 
5 FIGURE 12 are for illustrative purposes only. The user may select a number of other 
gradients to compare, and may select a variety of other sample intervals. 

The present invention also allows for 'cut' 1302, 'copy' 1304, 'paste' 1306, 'invert' 
1308, and 'flip' 1310 procedures to further expedite sequence design and modification as 
may be better understood with reference to FIGURES 1 3 A and 1 3B. FIGURE 1 3A simply 

10 illustrates the user's access to the abovementioned procedures by selection of the Edit menu 
item 350. FIGURE 13B provides an illustrative example of the various plots associated 
with a given sequence and an exemplary illustration of a copy procedure to expedite a 
sequence modification. On the PE plot 1320, a waveform has been selected for copying as 
indicated by a shaded selected area 1322. Generally, selecting an area is accomplished 

15 through clicking and dragging the desired area 1322 with the pointer device 140. The 
contents of the selected area 1322 are next copied to the computer 110 memory by 
choosing the Edit menu item 350 followed by the Copy 1 304 command from the Edit drop 
down menu 1312 of FIGURE 13A. Once the contents are copied to the computer 1 10 
memory, the user may then select a point in one of the displayed plots, for example, point 

20 1 324 at which to paste the copied waveform. As illustrated, the waveform pasted at point 
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1324 is an inversion of the original copied waveform. This is accomplished by selection of 
the invert command 1 308 from the Edit menu item 350 after the paste procedure. Thus, by 
this method, an FSE sequence, for example, can be designed from duplicating the parts of a 
generic SE sequence. It should be noted that, as illustrated, it is not necessary that a 
5 waveform be copied from the same plot. 

An important innovative aspect of the present invention is the capability for real- 
time communication between the MRI scan controller 130 and the scanning hardware, e.g., 
MRI unit 150, allowing for immediate design modifications and corresponding visual 
feedback. In the current preferred embodiment, seven binary files corresponding to the 

10 graphic waveforms on the screen have been generated. The seven binary files correspond 
to the preferred seven DACs, the gradient DACs 240, 242, and 244, the RF shaper DACs 
250 and 252, and the synthesizer DACs 231 and 232 interfacing the MRI control system, 
the computer 110, with the MRI scanning hardware. Thus, two binary files exist for 
synthesizers, two for shapers, and three for the gradients. Utilizing seven binary files for 

1 5 real-time conmiunication between the MRI hardware and control system is only a preferred 
embodiment of the present invention. As previously mentioned, the system and method are 
not, however, limited to such an arrangement but can be extended or reduced to any 
number of DACs depending upon the actual system on which the present invention is 
applied. 
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In order to provide real-time design and feedback, an accurate time frame reference 
is needed to be established. Since all of the digital electronic devices in the overall MRI 
system accept integer values generally limited by the bit-size of the associated DAC, small 
errors can accumulate during operation due to round-off of input data from various 
5 assignments and calculations when the sequence is built. This may sometimes lead to 
serious consequences for proper realization of a MRI pulse sequence. Thus, it is 
particularly important for a real-time interactive system, as described herein, to provide 
verification and round-off correction procedures between the displayed graphics and all the 
DACs to ensure time alignment throughout the entire sequence design process. The timing 

10 error due to roundoff depends on the settings of the time-resolution, generally on the order 
of microseconds, in the hardware. However, the accuracy of all calculations and roundoff 
subroutines in the user interface is preferably on the order of a single nanosecond. Every 
piecewise segment of the waveform is rounded off to be an integer-multiple of the three 
different time units, i.e., the unit for the RF pulses, the unit for the gradient waveforms and 

15 the unit for the sampling rate. Potential conflicts among the three time coordinates are 
resolved prior to the integral conditions, i.e. the gradient waveform specifications required 
to produce an echo, being applied. In a preferred embodiment of the present invention, 
satisfaction of the integral conditions are performed automatically by the underlying 
algorithm. This occurs not only at the initial sequence design stage, but also during 

20 gradient waveform modification, for example. Thus, the present invention provides a 

29 



• 



PATENT APPLICATION 
Docket No. 459713-1 

dynamic response to the user's adjustment by dynamic calculation and adjustment for 
satisfaction of the requisite integral conditions, e.g., by input via mouse 140 such as 
dragging of the gradient for modification thereto or by delta tuning. 

Satisfaction of the integral conditions, as had by the present invention, may better 
5 be understood with reference to FIGURES 14A and 14B. To obtain the largest echo 
amplitude, both the SS and RO gradients must satisfy the corresponding integral 
conditions. For the SS integral condition, the time integrated area of SS area 1400 must be 
equivalent to the time integrated area of SS area 1402. This is because the protons excited 
by the RF pulse 1404, dephase during the time period of SS area 1400 and subsequently 

10 rephase during the span of time corresponding to the SS area 1402. Similarly, the time 
integrated area of RO area 1406 must be equivalent to the RO area 1408. 

To refocus the spins that are in motion, higher order integral conditions are 
required. This scenario is illustrated in FIGURE 14B. In this illustration, a gradient 
sequence that has both stationary spins and spins moving in constant velocities are 

1 5 refocused at the time of the echo center, i.e. the signal components contributed by the spins 
in constant velocities are included in the image instead of, as would occur otherwise, the 
motion artifacts. Calculation of the integral condition for a particular sequence comprises 
the solution of a set of simultaneous integral equations depending on the order of 
compensation. As noted above, these integral conditions are calculated, and the sequences 
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are modified accordingly, both at the initial sequence design and dynamically as 
modifications are being made to a sequence. 

In FIGURE 15 A is illustrated a gradient sequence with velocity motion 
compensation for both the SS and RO gradients. The default timing settings of SS nodes 
5 1 500 and 1 502 and RO nodes 1 504 and 1 506, are also illustrated. For the SS gradient, the 
time and amplitude of node 1500 are initially 7.232 ms and -1 16.689, respectively, while 
the time and amplitude of node 1502 are initially 8.256 ms and 58.358, respectively. The 
SS gradient was then modified by dragging node 1502 to a timing position of 10.000 ms, 
denoted node 1502B in FIGURE 15B. In order to satisfy the integral condition such that 

10 the resulting echo will be formed at a substantially identical position as an echo resulting in 
the non-modified SS gradient as illustrated in FIGURE 15 A, the SS gradient must be 
modified. This modification is apparent in FIGURE 15B by the amplitudes of nodes 
1500B and 1 1 lOB automatically adjusting to -95.313 and 90.017, respectively so that the 
integral conditions are still satisfied. Similarly, in FIGURES 15A and 15B, node 1504 of 

15 the RO gradient has been dragged from an initial time of 6.720 ms to 5.536 ms. The 
immediate response displayed in the user interface is an adjustment of amplitude of RO 
node 1504 from 165.123 to 215.023 (indicated as node 1504B) and an adjustment of RO 
node 1506 amplitude of -279.913 to -241.386 (indicated as node 1506B). 



31 




PATENT APPLICATION 
Docket No. 459713-1 

In FIGURES 16A through 16D is illustrated a spin echo sequence example 
depicting motion compensation of the SS gradient for acceleration and RO for velocity 
compensation. 

As described, the present invention is generally composed of two portions: the 
5 graphical user-interface and the real-time MRI interface, or front-end. The real-time MRI 
interface allows communication between the graphical user-interface, and thereby the user, 
with the system hardware. The primary task of the front-end is to translate the output, e.g., 
the seven binary files, of the graphical user-interface into the data format required for input 
to the hardware and MRI controller and to retrieve the digitized signal from the system for 

1 0 either display in proper graphical form or delivery back to the user-interface for interaction 
purposes. Three types of analyzers for handling the three different kinds of signals, i.e. 
MRI signals, gradient waveforms and RF shapes, received from inside the aperture 155 are 
required. All three analyzers require the real-time feedback and capabilities for a fully 
real-time and interactive MRI process. 

15 Accordingly, a second window (in addition to the aforedescribed primary user 

interface 300 window, and the various manifestations thereof) is created when the real-time 
interface is invoked and is illustrated in FIGURE 17. 

Second window 1700 contains various menu editor items, e.g.. File 1702, Scan 
Control 1704, Seq Tuning 1706, Tuning Tools 1708 and Scan 1710. When a sequence is 

20 created along with the corresponding .exam and .va files, as described hereinabove, the 
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sequence can be selected by the options available (not shown) through the File menu editor 
item 1702. Li order to generate and receive the MRI signal, several system parameters 
have to be tuned properly for the loaded sequence. These parameters include, for example, 
the settings of the central frequency of the magnet 170, the power gain of the RF amplifier, 
5 the gain of the receiver, etc. Tuning of these parameters is performed through a scan 
control dialog box which is invoked by selection of the Scan Control menu editor item 
1704. The scan control dialog box preferably has numerous edit-boxes and selection 
buttons corresponding to the aforementioned parameters settings. Graphical displays of the 
signal intensify, spectrum and phase of the MRI echo are also preferably provided in 

10 subframes for easy user viewing. 

Once these parameters are tuned, calibration of the sequence may be initiated by 
user selection of the Seq Tuning menu editor item 1706 which preferably generates a 
sequence tuning dialog box 1800, including four control frames 1802, 1804, 1806, and 
1 808, for graphical display, as illustrated in FIGURE 1 8 A. Numerous controls are provided 

15 and are generally divided into a control menu section 1810 and control mode section 1812. 
Various radio buttons under the control menu section 1810 provide selection options for 
specific signals, e.g., input/output from the scanner, real and imaginary parts, and 
amplitude and phase of the MRI echo to be displayed in the upper and lower frames. The 
control mode section 1812 provides control to the user of, for example, sequence repetition 

20 times, toggling of each of the gradients, selection of the warp level of the PE gradient, etc. 
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In FIGURE 18B is illustrated an exp5iiplary artificial harmonic signal as may be 

displayed in frames 1802, 1804, 1806^d 1808 of the sequence tuning dialog box 1800. 

The real and imaginary part of the artificial harmonic are respectively plotted in frames 

1802 and 1806 while the M;^tude and phase of the echo are respectively displayed in 

5 frames 1 804 and 1 808. Jm addition to the MRI echoes, gradient waveforms and RF shapes 

O could also be disgj^^ed in the sequence tuning dialog box 1800 frames. 

i-D A scan parameters dialog box 1900, as illustrated in FIGUR^ 19, is generated by 

I n ^^^^^ 

5^' user selection of the scan menu editor item 1710, locatedjjrfTGURE 17, for performing a 

'''r^'. phantom scan as part of the calibration proces^^^^merous scan parameters are displayed 

10 therein, e.g., TE 1902, number of sjie^ 1904, TR 1906, NEX 1908, levels 1910, and slice 
plane 1914. Some value^^etfe modifiable by the user while static values are shaded so that 
the user can vipj^out not modify, the parameters. Selection of the scan button 1916 will 
initi^t^ascan, resulting in image reconstruction and display for user analysis. 



Preferably, there will be two modes of communication between the graphical user- 
15 interface and the front-end. A manual load is available any time a sequence is selected via 
the procedure described with reference to the second window 1700 of FIGURE 17. The 
sequence will be repeatedly transmitted to the MRI hardware and executed by the MRI 
scanner. Returning signals are accordingly displayed on sub-windows, or frames, in the 
Sequence Tuning dialog box 1800 of FIGURE 18. Accordingly, in this mode, the loaded 
20 sequence that is repeatedly sent to the MRI hardware remains constant during the scan, i.e. 
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the user is unable to modify the sequence once it is loaded and running, and thus the 
aforedescribed techniques for sequence modification are unavailable in this mode. 
However, changes made to the sequence can be loaded through the reload feature in the 
editor sequence tuning, as shown in FIGURE 17, particularly the Sequence Tuning menu 
5 editor item 1706. 

A second mode of communication between the graphical user-interface and the 
front-end is referred to as auto load and is initiated as soon as the calibration check box 610 
of the sequence tailor dialog box 600 is selected. Thereafter, any change of the sequence 
through the various sequence modification techniques are provided for the graphical user- 

10 interface and automatically reloaded and transmitted to the hardware. These resulting 
signals are immediately retrieved thereby providing the user with a real-time display of the 
various effects of the user modifications. 

Since the desired sequence often becomes distorted inside the aperture 155, the 
auto-load mode is particularly suitable for dynamic correction of the gradient waveforms. 

1 5 For instance, the phase information of an MRI echo may be calculated and provided to the 
graphical user-interface through the front-end. This information may be used to reset the 
reference frequency and phase of the sequence. Such an automated iterative procedure 
may be used for phase adjustment and alignment of the MRI echoes as well. 

Furthermore, initiation of the auto-load mode had by selection of the calibration 

20 check box 506 of the sequence tailor dialog box 600 preferably terminates the 
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aforedescribed automatic integral conditioning thereby providing more freedom for 
sequence manipulation. Concurrently, all related restrictions for modification are 
terminated. 

Although preferred embodiments of the method and apparatus of the present 
invention have been illustrated in the accompanying Drawings and described in the 
foregoing Detailed Description, it will be understood that the invention is not limited to the 
embodiment disclosed, but is capable of numerous rearrangements, modifications and 
substitutions without departing from the spirit and scope of the invention as set forth and 
defined by the following claims. 
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